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' We present a microscopic theory of coherent quantum transport through a superconducting film 

, between two ferromagnetic electrodes. The scattering problem is solved for the general case of 

' ferromagnet/superconductor/ferromagnet (FSF) double-barrier junction, including the interface 

transparency from metallic to tunnel limit, and the Fermi velocity mismatch. Charge and spin 
conductance spectra of FSF junctions are calculated for parallel (P) and antiparallel (AP) align- 
ment of the electrode magnetization. Limiting cases of nonmagnetic normal-metal electrodes (NSN) 
' and of incoherent transport are also presented. We focus on two characteristic features of finite 

size and coherency: subgap tunneling of electrons, and oscillations of the differential conductance. 
Periodic vanishing of the Andreev reflection at the energies of geometrical resonances above the 
superconducting gap is a striking consequence of the quasiparticle interference. Also, the non-trivial 
spin-polarization of the current is found for FSF junctions in AP alignment. This is in contrast 
with the incoherent transport, where the unpolarized current is accompanied by excess spin accu- 
mulation and destruction of superconductivity. Application to spectroscopic measurements of the 
^ ' superconducting gap and the Fermi velocity is also discussed. 
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I. INTRODUCTION 



During the past decade, there has been a growing interest in various electronic systems driven out of equilibrium 
by the injection of spin-polarized carriers. Such systems can be realized by current-biasing structures consisting of 
ferromagnetic and non-ferromagnetic (e.g. superconducting) layers, due to the difference in population of majority 
and minority spin subbands.EI The concept of spin-polarized current nowadays has attracted considerable interest in 
^ ' ferromagnetic heterostructures, in particular for applications in spintronics.B 

Charge transport through a normal metal/superconductor (NS) junction, with an insulating barrier of arbitrary 
strength at the interface, has been studied by Blonder, Tinkham, and Klarwijk (BTK),0 and the Andreev reflec- 
tion is recognized as the mechanism of normal-to-supercurrent conversion. &B The BTK theory has been extended 
Xf^ by Tanaka and Kashiwaya to include the anisotropy of the pair potential in d-wave superconductors.ljLl The mod- 
t ification of the Andreev reflection by the spin injection from a ferromagnetic metal into a superconductor in fer- 
i romagnet/superconductor (FS) junctions was first analyzed by de Jong and Beenakker.l3 More recently, the effects 
' of unconventional d-wave and p-wave pairing and of the exchange interaction in FS systems_such as the zero-b 
, conductance peak and the virtual Andreev reflection, have been clarified by Kashiwaya et alu and Yoshida et al 
The Fermi velocity roismatch between two metals can also significantly affect the Andreev reflection by altering the 
-j->; , subgap conductance,Ej which is similar to the presence of an insulating barrier £3 

' In experiments, a superconductor is used to determine the spin polarization of the current injected from (or into) a 
ferromagnet by measuring thej-diffierential conductance. Thpsj? measurements have been performed on tunnelimctions 
in an external magnetic fieldjljo metallic point contacts, E£lE3 nano-contactaJormed by microlithography,EZl and FS 
■ O junctions with d-wave superconductors, grown by molecular beam epitaxy.ll3_[n diffusive FS junctions, the_epa;ess 
5 resistance may be induced by spin accumulation near the insulating interfacejlj and by the proximity effect. EjOcj 
When electrons pass incoherently through the interfaces, the BTK model can be successfully applied to 
normal naetal/superconductor /normal metal (NSN) or ferromagnet/superconductor/ferromagnet (FSF) double 
junctions.EScJ However, the properties of coherent quantum transport in clean superconducting heterostructures are 
strongly influenced by size effecjtST-|Which are not included in the BT K . ni o d el. Well-known examples are the current- 
carrying Andreev bound stateaD'Cj and multiple Andreev reflectkmsr^rlrl in superconductor/normal metal/super- 
d ' conductor (SNS) junctions. Since early experiments by Tomasch,t3 the geometric resonance nature of the differential 
conductancej QSC |i ll ajt4pns in SNS and NSN tunnel junctions has been ascribed to the electron interference in the 
central film.E3'E2ll23'E2l Recently, the McMillan-Rowell oscillations were observed in SNS edge jjiiictions of d-wave 
superconductors, and used for measurements of the superconducting gap and the Fermi velocityEj 
Here we presenL.a_comperhensive microscopic theory of coherent transport in FSF double junctions (with NSN as 
a special case)E3'L3 We limit ourselves to . clea n conventional (s-wave) superconductors, and neglect, for simplicity, 
the self-consistency of the pair potential,EZI'E3 and nonequlibrium effects of charge and spin accumulation at the 



* E-mail: mbozovicainfosky.net 
t E-mail: zradovic9ff.bg.ac.yu 



2 



interfacesrjo When two interfaces are recognized by electrons simultaneously, characteristic features of finite size 
and coherency are the subgap tunneling of electrons and oscillations of both charge and spin differential conductances 
above the gap. One consequence of the quasiparticle interference is the periodic vanishing of the Andreev reflection 
at the energies of geometrical resonances. The other is the existence of a non-trivial spin-polarization of the current 
not only for the parallel (P), but also for the antiparallel (AP) alignment of the electrode magnetizations. Previous 
analysis of incoherent transport in FSF double junctions predict the absence of the spin current and suppression of 
supe g^^ ductivity with increasing voltage for AP alignment, as a result of spin imbalance in the superconducting 



II. THE SCATTERING PROBLEM 



We consider an FSF double junction consisting of a clean superconducting layer of thickness /, connected to fer- 
romagnetic electrodes by thin, insulating interfaces. Fig. 1. For the ferromagnetic metal we adopt the Stoner 
model describing the spin-polarization effect by the usual one-electron Hamiltonian with an exchange potential. The 
quasiparticle propagation is described by the Bogoliubov-de Gennes equation 



i7o(r)-p./i(r) A(r) 



*,(r) = ^^^(r), (1) 



with Hq{v) = -h^V"^ /2m+W {y)^U {y)- ^i, where C/(r) and /i are the Hartree and the chemical potential, respectively. 
The interface potential is modeled by W{y) = W{5{z) + S(z — I)}, where z-axis is perpendicular to the layers and 
6{z) is the Dirac J-function. Neglecting the self-consistency of the superconducting pair potential, A(r) is taken 
in the form A<d{z)<d{l — z), where A is the bulk superconducting gap and 0(z) is the Heaviside step function. In 
Eq. (|^), cr is the quasiparticle spin {a =t, | and a =|, t), E is the energy with respect to ^, h{r) is the exchange 
potential given by ho{Q{—z) + ]8(2: — I)} for the P [AP] alignment, and pa is 1 (-1) for spins up (down). The 
electron effective mass m is assumed to be the same for the whole junction. Here, fj, — U{r) is the Fermi energy of 

the superconductor, E^p\ or the mean Fermi energy of a ferromagnet, E^p"^ = {Ep + Ep)/2. Moduli of the Fermi 



wave vectors, fc^^'' = \J 2mE'^p^ / and k!^p^ = y2mE^p'' /h?, can be different in general, and in the following, the 

Fermi wave vector mismatch (FWVM) will be taken into account through the parameter n — fcj^-* / k^p^ . The parallel 
component of the wave vector k|| a- is conserved, and the wave function 

*^(r) = exp(ik||^<, • r) ^^(z), (2) 



satisfies the boundary conditions 



tpa{z)\z=0- = 1pa{z)\z=0+, (3) 



dil)a{z) _ dil)c{z) 



dz 



dz 
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-^^^(0), (4) 
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FIG. 1: Double barrier junction consisting of two ferromagnets (F) and a superconductor (S) separated by insulating bar- 
riers (I), for parallel (left panel) and antiparallel (right panel) alignment of the electrode magnetization. Schematic of the 
corresponding densities of states. 
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1p,Tiz)\z=l- = i'criz)\z=l+, (5) 



dz 



dz 



-jr^Ml)- (6) 



Four independent solutions of Ea. (|l|) correspond to the four types of injection: an electron or a hole from either the 
left or from the right electrode.cl 

For the injection of an electron from the left, with energy E > 0, spin a, and angle of incidence 6 (measured from 
the z-axis), solution for ^a{z) in various regions has the following form: 
in the left ferromagnet (z < 0) 

4^,{z) = {exp(zA:+z) + b,{E, 0) exp(-ifc+z)} (^j^^+a^{E, 0) exp{ik^ z) J j , (7) 

in the superconductor {0 < z < I) 

Mz) = {ci(£;,0)exp(jg+z) + C2(£;,0)expHg+z)} (^Jf 



+{C3{E, 0) exp{iq^ z) + c^{E, 0) exp(-ig^ z)} [Jy* ) ' 
and in the right ferromagnet (z > Z), for the P [AP] alignment of the magnetizations 

Mz) = cAE, 0) eMKl.f) [l)+da{E, 0) eM-^Kl^f) ( 1 ) " (9) 

Here, u ^ + n/E)/2 and v = - n/E)/2 are the BCS coherence factors, and n = ^fW~~^ . The z- 
components of the wave vectors are fc^ = [2'm/h?){E^p^ + pah^ zL E) — k^| ^, and = 

y^(2m/ft2)(i;(_f) ±n)- k^^^, where |k||,^| = {2m/h?){E^/^ + p^iQ + E) sm0. The coefficients a^, and 

c?cr are, respectively, the probability amplitudes of: (1) Andreev reflection as a hole of the opposite spin (AR); (2) 
normal reflection as an electron (NR); (3) transmission to the right electrode as an electron (TE); (4) transmission 
to the right electrode as a hole of the opposite spin (TH). Processes (1) and (4) are equivalent to the formation of a 
Cooper pair in the superconductor by taking one more electron from either the left or the right electrode, respectively. 
Amplitudes of the Bogoliubov electron-like and hole-like quasiparticlcs, propagating in the superconducting layer, 
are given by the coefficients ci through C4. 

From the probability current conservation, the probabilities of outgoing particles satisfy the normalization condition 

AAE,0) + B^E.ff) + CAE,0) + DAE,0) = 1, (10) 

where, 

A^E^O) = K{EM\ (11) 

BAE,0) = K{E,e)\\ (12) 
C^E^O) = "^(^f^ \ca{E,6)\\ (13) 

D,{E,6) - ^(h^\d,[E,0)\\ (14) 

Neglecting small terms E/E^j^^ <C 1 and A/i?jf'' ^ 1 in the wave vectors, except in the exponents 

C±^l{q+±q-), (15) 
solutions of Eqs. (||)-(|^) for the probability amplitudes can be written in the following form 

aAE,0) = 4(fc./g.)Asin(C-/2) [^^fl j^,i^(^_/2) + zg^l] cos(C-/2)] , (16) 
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b,{E, 9) = -[A^C+A^ - {A^C+E^ + B^V+Q^) cos(C-) + {A'^C- + B^V^) fl^ cos(C+) 



{B^C++A^V+) Ensm{C_) - i {B^C- + A":V^) sin(C+)], 



(17) 



d,{E, 



where 



r 

x{i [J^+ cos(C+/2) + i£+ sin(C+/2)] £;sin(C_/2) - [£+ cos(C+/2) + iJ^+ sin(C+/2)] f7 cos(C-/2)}, (18) 
4(fc^/gJAr!e'^^' ^ 

r ^ 

XI [T^ cos(C+/2) + sin(C+/2)] sin(C-/2), 

r = A^A'^A'^ - (AiA'^E^ + B^B'^n^) cos(C-) + {A^A^ + B^B^) cos(C+) 
+i {A^B'l + B^A^) EilsiniC^) ~ I {A^_B^1 + B^_A^_) V? sin(C+). 



(19) 



(20) 



In Eqs. ([|)-(|^) 



A 



L{R) 



L(R) 



B 



L{R) 



1 ± i^f 



C±=K^ TK^: 

^± = 1 ± K^K^, 



with K 



V± = -(1 T K^*Kl^), £± = ± Ki, 

{ = (k^ + /q,, - (k^ - iZ) /q,, Kf = (fc,[^] + zz) /q,, = (fc^i^j - /g„ for the P [AP] 



ahgnment. Here, K^* — {kc — iZ)/qa is the complex conjugate of , and Z ~ 2niW /h?k^p^ is the parameter 
measuring the strength of each interface barrier. Approximated wave- vector components, in units of kp \ are 



(ja = yl ^ ~ Ao-cosS, and |k||,cr| — AcrSin^, where Ao- — K^/l + paX, X = ho/E)^' > 0, and k 7^ 1 is 

measuring FWVM. 

In the corresponding FNF double junction, AR and TH processes are absent, and the expression for NR amplitude, 
Eq. (p^), reduces to 



(Kt ~ Kj^) cos(/g^) + z(l - Kt Kj") sin(;g^) 
[K^ + K^) cos{lq^) - z(l + K[K^) sm{lqN) ' 



(21) 



N 



(2m/ft2)(45)+i?)-k2^^. 



where 

To complete our considerations, we also present the probability amplitudes for an PS single junction in the same 
notation, 



da{E,( 



2(fc./g.)A 
A^E + B^n' 
C+E + v+n 
A^E + B'^n' 

2{K/q,)Eu{\ + K^ 

A^E + B^Q. 
2{Klq,)Ev{\ - Kj 

AiE + Bin 



(22) 
(23) 
(24) 
(25) 



Note that Co- and now describe the transmission of the Bogoliubov electron-like and hole-like quasiparticle, xe- 
spectively. The well-known BTK results can be reproduced by taking X = 0, k = 1, and = in Eqs. (2^-(^).Ej 
In case of the simplest NSN metallic junction^^aking X ^ 0, Z = 0, and k = 1 in Eqs. (p^-(p^), the scattering 
probabilities can be written in an explicit formE3 



Aa{E,9) = 

Ba{E,e) = 

Ca{E,9) - 



Asin(C-/2) 



£;sin(C_/2) + irJcos(C-/2) 
D„{E,9)^0, 



£;sin(C-/2) + ir2cos(C-/2) 



(26) 
(27) 
(28) 
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Solutions for the other three types of injectfon can be obtained by the same procedure. In particular, if a hole with 
energy —E, spin cr, and angle of incidence 9 is injected from the left, the substitution ^ q~ holds, and the 
scattering probabilities are the same as for the injection of an electron with E, a, and 9. Therefore, in order to 
include the description of both electron and hole injection, the calculated probabilities should be regarded as even 
functions of E. Also, for an electron or a hole, injected from the right, the probabilities are the same as for the 
injection from the left, except a —^ a for the AP alignment. 

Following the conservation of k|| o-, transmission of an electron (hole) with a =t, injected from the left electrode into 
the superconductor, is possible only for angles of incidence 9 satisfying 9 < 9ci, where 9ci = arcsin(l/Af ) is the angle 
of total reflection. Then, Ai{E,9) ~ and B^{E,9) = 1 for 6* > 9ci- On the other hand, ki, which corresponds to 
the hole (electron) created by the Andreev reflection, is real only for 9 < 9c2 — arcsin(A|/A|). The virtual Andreev 
reflection occurs for 9c2 < 9 < 9ci, since becomes imaginary in that case. For injection of an electron (hole) with 
a =1, transmission into the superconductor is possible for any 9 < tt/2, and fc| is always realJ3 
From Eqs. (H) and (|l|) it follows that A^{E, 9) ^ D„{E, 9) ^ Q when 

C_ = 2mi (29) 

for n — 0, ±1,±2, . . .. Therefore, the Andreev reflection at both interfaces vanishes at the energies of geometrical 
resonances in quasiparticle spectrum. The effect is similat-to the over-the-barrier resonances in the simple problem of 
one-particle scattering against a step-function potential,cHI the superconducting gap playing the role of a finite-width 
barrier (as in the semiconductor modelcj). The absence of AR and TH processes means that all quasiparticles with 
energies satisfying Eq. (^9|) will pass unaffected from one electrode to another, without creation or annihilation of 
Cooper pairs. 

Characteristic features of coherent quantum transport through clean superconducting layers are the subgap tunneling 
and oscillations of the scattering probabilities. For E < A, the subgap tunneling suppresses the Andreev reflection, 
thereby enhancing the transmission. For E > A, all probabilities oscillate with E and / due to the interference effect. 
The interface resistance reduces AR and TE, and enhances NR and TH probabilities. In contrast to the positions of 
zeros of A^{E,0), given by Eq. (p9|), the positions of maxima of Ao-(£',0), as well as that of zeros and maxima of 
Bcr{E, 0), Ca-{E, 0), and D„{E, 0), are Z-dependent. Approaching the tunnel limit {Z oo), peaks in the scattering 
probabilities gradually split into two spikes belonging to consecutive pairs with positions defined by the quantization 
conditions 

= niTT, lq~ = n27r. (30) 

Here, ni — n2 — 2n, with n coresponding to that of Eq. (p9|). The exception is the spike at the gap edge, originating 
from the singularity in the BCS density of states. Note that Eq. ( ^0|) gives the bound state energies of an isolated 
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FIG. 2: Scattering probabilities, Aa{E,0), Ba{E,0), Ca{E,0), and Da{E,0), for an NSN double junction with transparent 
and highly non-transparent interfaces, Z = (dotted curves) and Z — 10 (solid curves), respectively. The parameters are 
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FIG. 3: Scattering probabilities, Ac{E,0), Bcr{E,0), C'cr{E,0), and Dcr{E,0), for an FSF double junction with thin supercon- 
ducting film, Ik^p^ = 10^, for X — 0.5, Z = 0, k = 1, A/E^p^ — 10~^, and P alignment. Solid curves: injection of an electron 
with a =t. Dashed curves: injection of an electron with a =|. 
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FIG. 4: Scattering probabilities, Aa{E, 0), Ba{E, 0), Ca-{E, 0), and Da{E, 0), for an FSF double junction with thick supercon- 
ducting film, Zfc^' = 10'*, for X = 0.5, Z = 0, k = 1, A/E^p^ — 10~^, and P alignment. Solid curves: injection of an electron 
with a =t. Dashed curves: injection of an electron with a =j. 



superconducting film. This is illustrated in Fig. ^ for an NSN junction, showing a simple connection between the 
resonances in metallic junctions {Z — 0) and the bound states in the corresponding tunnel junctions (Z — > oo). 
The influence of exchange interaction is illustrated in Figs. |^ and ^ for an FSF double junction in P alignment with 
Z = and K = 1. Taking A/£'jf^ = lO'^, in a thin superconducting fllm, /fc^-* ~ lO'^, the Andreev reflection 
is strongly suppressed, since the subgap transmission is considerable. Fig. 0. In this case, the oscillations are less 
pronounced, with the period much larger than A. For a thick film, Ikp ^ 10 , the subgap tunneling is irrelevant 
(except for small 'tails' in Aa-{E, 0) and Ccr{E, 0) at E < A) and above the gap the oscillations arc pronounced, with 
the period on the order of A, Fig. ^. The scattering probabilities for AP alignment differ very slightly in the case 
of normal incidence, 6 = 0. Although spin-independent due to the singlet-state pairing, Aa{E,0) is suppressed in 
comparison with the corresponding NSN junction, and D„{E,0) becomes non-trivial. The spin-dependent normal 
reflection also occurs, B„{E,Q) having zeros at the same energies as 0) and Dcr{E,0), so that maxima in 

C'a-{E,0) are still equal to unity due to the interface transparency. 
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III. DIFFERENTIAL CONDUCTANCES 



When voltage V is applied to the junction, the charge current density is given by 

h{y) = Y.j (^^^•^'^/(k-^): (31) 

where v = (?i/m)S[u* (r)VMcr(r) + uJ.(r)Vwo-(r)] is the velocity, and (5/(k, V) is the asymmetric part of the nonequi- 
librium distribution function of current carriers. Using the solution of the scattering problem for the injection of an 
electron from the left, described in the previous section, and the dispersion relation k(i?) = A;+z + k|| o-, Eq. ( |3l] ) can 
be rewritten in the form 

hi^)^^^ / dEY,K / dBsmBcosB[l + A,{E,e) - B,{E,e)]5f{\<i,V). (32) 

-cx) 

In accordance with BTK, without solving the suitable transport equation, we take (5/(k, V) = fa{E — eV/2) — fo{E + 
eV/2), where fo{E) is the Fermi-Dirac equilibrium distribution function. 
In this approach, the charge current per electron is given by 

oo 

I,iV)^^ J dE[fo{E~eV/2)- fo{E + eV/2)]Gg{E), (33) 

— oc 

where the spin-averaged differential charge conductance at zero temperature is 

n/2 

Gg {E) = —J2>^1 / dd sin 9cos9[l + A^{E,e)-B,{E,d)]. (34) 




FIG. 5: (Left panel) Differential charge (a) and spin (b) conductance spectra, Gq{E) and Gs{E), of an FSF double planar 
junction with thin superconducting film, Ik^p^ = 10^, for X — 0.5, Z = 0, k — 1, A/E^p^ = 10~^, in P and AP alignment. 
Conductances of the corresponding NSN junction are shown for comparison (dotted curves). Arrows indicate and G^ 
values. 

FIG. 6: (Right panel) Differential charge (a) and spin (b) conductance spectra, Gq{E) and G3{E), of an FSF double tunnel 
junction with thin superconducting film, Ik^p^ — 10'', for X — 0.5, Z = 10, k = 1, A/E^p'' — 10~^ , in P and AP alignment. 
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FIG. 7; (Left panel) Differential charge (a) and spin (b) conductance spectra, Gq{E) and Gs{E), of an FSF double planar 
junction with thick superconducting film, Ik^p^ = 10^, for X — 0.5, Z = 0, n = 1, A/E^p '^ = 10""^, in P and AP alignment. 
Dashed curves represent the generalized BTK results for the same parameters. Arrows indicate and G^ values. 

FIG. 8: (Right panel) Differential charge (a) and spin (b) conductance spectra, Gq{E) and Gs{E), of an FSF double planar 
junction with thick superconducting film, Ik^p^ = 10*, for X = 0.5, Z = 1, k = 1, A/E^p^ = 10"'^, in P and AP alignment. 
Dashed curves represent the generalized BTK results for the same parameters. Arrows indicate Gq and values. 



By analogy, the corresponding spin current (proportional to the probabihty current) is given by 

oo 

Is{V) = \ j dE [fo{E - eV/2) - fo{E + eV/2)] Gs{E), (35) 



where the differential spin conductance at zero temperature is 

t/2 



2 /> 

GsiE)^ ^Y^p^Xl d0 sm9cos9[l- A,{E, 9)- B,iE,e)]. (36) 



Note that in Eqs. (|3J) and ( |36| ) the upper limit of integration over 9 for a =1 is the angle of total reflection 9ci < tt/2. 
Avoiding integration over 0, the differential charge and spin conductances of a point-contact FSF double junction 
are simply expressed by 



and 



^^(^) = I^E^-^' [1 - ^a(£;,0) - B^iE,0)] . (38) 



The influence of the exchange interaction on the conductance spectra is shown for X — 0.5 on the example of thin, 
= 103, and thick, Ik^^ = lO'', superconducting films (Figs. |[^). Besides the case of transparent interfaces, 
Z = (Figs. 1^ and the effect of interface resistance is illustrated in the tunnel limit, Z = 10, and for weak 
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FIG. 9: (Left panel) Differential charge (a) and spin (b) conductance spectra, Gq{E) and Ga{E), of a point-contact FSF 
double junction, with P alignment, for the same parameters as in Figs. 7 and 8 {X = 0.5, A/E^p^ = 10~^, Ik^p ^ = 10'*, and 
K = 1). Here, the difference between conductances for the AP and the P alignment is negligible. 

FIG. 10: (Right panel) Square of resonant energies E^ as a function of n^, obtained from positions of minima of G{E) for 
Ik^p^ = 10''. The intercept is and the slope is {hv'-p^ /2lf . 

non-transparency, Z — 1, in FigSj. ^ and The influence of FWVM on the conductance spectra, k ^ 1, is similar to 
that of the interface resistance.Ea The values of normal conductances, and of the corresponding FNF double 

planar junction, indicated by arrows, are obtained by setting A^{E,e) = and B„{E,e) = \b^{E,e)\'^ in Eqs. (||) 
and (11), where {E, 9) is given by Eq. (|l]). 

The spin-polarized subgap tunneling of quasiparticles, and strong suppression of the Andreev reflection as a con- 
sequence, is significant in thin superconducting films, whereas the conductance oscillations above the gap are pro- 
nounced in the thick films. The magnetoresistance is apparent, as charge and spin conductances are larger for the 
P than for the AP alignment. An important consequence of the coherent transport is that the spin conductance 
is non-trivial for the AP alignment, approaching its normal value G^ — either for E/A. 3> 1, or in the tunnel 
limit [Z oo) for all energies. We emphasize that the amplitudes of the oscillations are considerably larger for the 
point-contact FSF than for the planar FSF double junction, Fig. ||. 

Incoherent transport through an FSF double junction is described as a transport through the corresponding FS and 
SF junctions in series. In that case, the conductance spectra are calculated using the generalized BTK probabilities, 
obtained from Eqs. ( [2^ ) and (p3[). Numerical results for the incoherent transport are also presented in Figs. |^and ||for 
comparison. It is evident that in thick films the only difference comes from the interference-effect oscillations for the 
energies above the gap. In contrast with the coherent transport, for the AP alignment Gs{E) = 0, and nonequilibrium 
spin density accumulation changes the chemical potential of two spin subbands in the superconductor. This reduces 
the supecconducting gap with increasing voltage, and destroys the superconductivity at a critical voltage on the order 
of A/e.r'^l.The effect of incoherency is less pronounced in metallic than in the tunnel junctions due to the Andreev 
reflection. Ej 

(s) 

The results can be applied to reliable spectroscopic measurements of A and Vp in superconducting films. From Eq. 
(p9|), for 9 — 0, the energy En of the n-th geometrical resonance (conductance minimum) satisfies a simple relation 




Therefore, the linear plot of E^ vs has the intercept equal to and the slope equal to {hvp /2l)'^. An example 
is shown in Fig. |l^. Note that even the points obtained for planar (3D) double junctions lie almost on the same 
straight line given by Eq. ( ^9| ) for the particular case of a point-contact (ID) double junction. The numerical results 
show that the method is almost independent on dimensionality of the junction and on parameters of the electrodes. 
The net spin polarization of the current is defined as n(V^) = Is{V)/ Iq{V). In thin superconducting films, n(F) is 
almost constant, considerably smaller than X, which is the polarization in the corresponding FNF junction. Below 
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the gap, ey/2A < 1, the contribution of subgap tunnehng to spin polarization becomes neghgibly smah for large 
I. Above the gap, eV/2A > 1, the polarization increases with V, the increase becoming steeper with the interface 
non-transparency. Fig. On the other hand, in a tunnel FS junction the polarization abruptly changes from 11 = 
to n = X at eV/2A = 1. The same result holds for incoherent transport through an FSF double tunnel junction. 
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FIG. 11: Spin polarization of the current, HCl/), for an FSF double planar junction with P and AP alignment, is shown for 
(a) Z = 0, and (b) Z = 1. Other parameters are the same as in Figs. | and | (X = 0.5, A/e'^^ = 10"^, Zfc^f' = lO"*, and 
K — 1). Dashed curves represent the generalized BTK results. 



IV. SUMMARY 

We have analyzed transport properties of FSF double-barrier junctions, taking into account the influence of the 
exchange interaction, the resistance of the interfaces, and the Fermi velocity mismatch on the scattering probabilities 
and the conductance spectra. We have shown that subgap tunneling and oscillations of differential conductances are 
the main features of the coherent quantum transport through a superconducting layer in both FSF and NSN double- 
barrier junctions. The subgap tunneling suppresses the Andreev reflection, thereby enhancing the transmission, 
especially in thin films. The scattering probabilities and conductances oscillate as a function of the layer thickness 
and of the quasiparticle energy above the gap. 

Periodic vanishing of the Andreev reflection at the energies of geometrical resonances is found as an important 
consequence of the quasiparticle interference. Insulating barriers at the interfaces reduce the Andreev reflection 
and transmission, mainly for energies below the gap. The Fermi velocity mismatch has a similar effect. Results are 
directly accessible to experiments. In principle, oscillations of differential conductances with the period of geometrical 
resonances could be used for reliable spectroscopy of quasiparticle excitations in superconductors. 
In conclusion, finite-size effects, along with the difference between coherent and incoherent transport, are essential 
for spin-currents in FSF junctions. For the coherent transport, besides the spin-polarized subgap tunneling in 
thin superconducting films, pronounced oscillations of spin conductance occur in thick films. As a consequence of 
the quasiparticle interference, a non-trivial spin current without the excess spin accumulation and destruction of 
superconductivity is possible even for AP alignment of the electrode magnetizations. 
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